Global dynamics of cosmological scalar fields - Part II 
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Abstract 

This is the second part of integrability analysis of cosmological models with scalar fields. 
Here, we study systems with conformal coupling, and show that apart from four cases, where 
explicit first integrals are known, the generic system is not integrable. We also comment on 
some chaotic properties of the system, and the issues of integrability restricted to the real 
domain. 

1 Introduction 

Conformally coupled fields were subject to more rigorous integrability analysis, as opposed 
to minimally coupled ones, thanks to the natural form of their Hamiltonian. As will be shown 
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in the next section, the kinetic part is of natural form, albeit indefinite, and the potential is 
polynomial (in the case of real fields). 

Ziglin proved that the system given by (|8j) is not meromorphically integrable when A = A = 
and k = 1 [18] . His methods were then used by Yoshida to homogeneous potentials which is 
the case for the system when k = [141 fl"5j [T6j 117] . Later, Yoshida's results were sharpened 
by Morales-Ruiz and Ramis [llj, and used by the present authors in [9] to obtain countable 
families (restrictions on A and A) of possibly integrable cases. Also recently, more conditions for 
integrability have been given in [2], although only for non-zero spatial curvature k and a generic 
value of energy, that is, when the particular solution is a non-degenerate elliptic function (in 
particular not for zero energy). 

Our work shows, that the conjecture of that paper is in fact correct - as shown in Section 4 
- the system is only integrable in two cases (with the above assumptions). We go further than 
that and show that for a generic energy value, a spatially flat (k = 0) universe is only integrable 
in four cases. Also, the particular case of zero energy is analysed and new, simple conditions on 
the model parameters are found. Finally, we check that when E = k = 0, the problem remains 
open, as the necessary conditions for integrability are fulfilled. 

When it comes to numerical studies of the problem, there are various results, most notably 
chaotic behaviour [5], but also fractal structure and chaotic scattering [13] . However, it remains 
unclear whether the widely exercised complex rotation of the variables changes the system's 
integrability. Even for very simple systems it was shown [3] that there might exist smooth 
integrals, which are not real-analytic. This question is especially vital since our Universe clearly 
has real initial conditions and dynamics. 

In what follows, we derive the Hamiltonian system for the conformally coupled scalar field, 
and proceed to analyse its integrability. For the basics of the theory used, and relevant, more 
detailed literature see the first part [T2] . 



2 Conformally coupled scalar fields 



The procedure of obtaining the Hamiltonian is the same as in the case of minimally coupled 
fields, only this time the action is 

^d 4 x, (1) 
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where an additional coupling to gravity through the Ricci scalar 1Z, and a quartic potential term 
with constat A are present, as opposed to the minimal scenario. The cosmological constant A 
and the mass m remain as previously. We keep the same notation as before and express the 
involved quantities in comoving coordinates and conformal time to get 



C = 6(a"a + Ka 2 ) - ^a"a\^\ 2 + ^'|V 



m 



A a >| 4 - 2Aa 4 - l -Ka 2 



(2) 



from which we remove a full derivative, and introduce new field variables ifi = yl2 (f>exp(i0)/a 
to obtain 



£ = 6 
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(3) 
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The associated Hamiltonian is 



2,1 2 



+ 6 



(4) 



We can see that 9 is a cyclical variable because we took the potential to depend on the modulus 
of ip only, so we write a constant instead of the respective momentum pg = uj. 

Finally, we express everything in dimensionless quantities, rescaling the constants, but also 
the time and momenta (as they are in fact time derivatives), which results in rescaling the whole 
Hamiltonian. We do this as follows 
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when K ^ 0, and using another of the dimensional constants otherwise. Thus, eliminating a 
multiplicative constant, the Hamiltonian reads 

,2 



H 



1 



{pI-pI) + 



1 



+ 



i 2 2 i.2 

+ mam 



+ - (Aa 4 + A0 4 ) 



(6) 



with k £ {-1, 0, 1} (K = k\K\)\ uj, A, A, m 2 £ R, and H = in any physically possible setup. 
However, the addition of radiation component, which scales like a -4 in the original action, leads 
to a constant in the Hamiltonian. Thus, it justifies studying energy levels other than zero as 
well. 

We note that for m = the system decouples, and is trivially integrable as shown in Appendix 
A. That is why we will assume m ^ henceforth. We will also take u = 0, that is, consider a 
scalar field equivalent to a real field after a unitary rotation in the complex ip plane. 

We change the field variables into the standard q and p ones for further computation, taking 



qi, 

12, 



Pa 

P<t> 



Pi, 
P2- 
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The Hamiltonian is then 

H 



^{-pI + pI)+v, 

V=\ [H-qj + q 2 2 ) + m 2 q 2 q 2 ] + - {Aqf + Ag 4 ) . 



(8) 



3 Known integrable families 

There are four known cases when the system has an additional first integral, independent 
of the Hamiltonian. They were found by applying the so called ARS algorithm basing on the 
Painleve analysis The following table summarises those results. 



solvability case 


k 


A 


2 

m 


(1) 


0,±1 


A = A 


m 2 = -3A 


(2) 


0,±1 


A = A 


m' 2 = —A 


(3) 





A = 16A 


m' 2 = — 6A 


(4) 





A = 8A 


m 2 = -3A 
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And the respective integrals of the systems are 
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In this work, we will show, that the above are the only integrable cases, when m 7^ 0. An 
important point to note is that there is a complete symmetry with respect to interchanging A 
and A. It is a consequence of the fact, that there exists a canonical transformation of the form 



Pi^ipi, qi^-iqi, 

P2^P2, 92^92, 



that changes the Hamiltonian to 

H = \ (Pi + Pi) + \ [k(qf +ql)- m\h 2 2] + \ (Mf + Agf) , 



(10) 
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which is the same after swapping the indices. We shall use this form of H, where the kinetic part 
is in the natural form, to make the use of some already existing theorems more straightforward. 



4 Integrability of the reduced problem 

It is possible to give stringent conditions for integrability of the system, by considering a 
reduced Hamiltonian. Namely, we can separate potential V into homogeneous parts of degree 2 
and 4: 

V = V h2 + V h4 , 

V h 2 = ^ k (9i + 9!) , (12) 

V M = - {-2m 2 q\q\ + Aqj + Xq 2 ) . 

The following fact is crucial in our considerations: if a potential V is integrable then its higheest 
order as well as the lowest order parts are also integrable. For the proof, see [5) This means that 
in our case if V given by (fT2j) is integrable then Vh 2 and must also be integrable. Vh 2 is the 
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potential of the two-dimensional harmonic oscillator, thus, it is trivially integrable. However, 
the homogeneous part Vh4 gives strong integrability restrictions for the whole potential V. We 
will call Vm the reduced potential and denote it by V. 

Thus we effectively set k = 0, and are now in position to exercise known theorems concerning 
homogeneous potentials of two variables. In particular the complete analysis for degree 4 has 
been completed in |10j . 

In order to identify our potential with some of the list given in that paper, we have to check 
how many Darboux points there exist, and what are the values of parameters A, A and m that 
give potentials equivalent to a particular family. 

We say that a non-zero point (qi, (72) = d is a Darboux point of the potential V(qi, (72) when 
it satisfies the equation 

V\d) = 7 d, (13) 
where 7 G C* = C \ {0}. Such a point corresponds to a particular solution of the form 

q(v) = f(v)d, p(v) = f(v)d, (14) 
with f{rf) satisfying the differential equation (for degree 4 potential) 

m = -ifivf- (15) 

As explained in the first part of this paper, particular solutions allow for studying the vari- 
ational equations along them, and yield necessary conditions for existence of additional first 
integrals. However, the major simplification discovered in [10J is that additionally there is only 
a finite number or parameter sets (or non-equivalent potentials) corresponding to integrable 
cases. 

Following the cited paper's exposition (and notation) we find that our potential has: 

1. Four simple Darboux points, when A(m 2 + A)(m 2 + A) 7^ 0, and AA 7^ m 4 . The only 
integrable cases are: 

(a) A = A = —\m 2 (V is equivalent to V4), 

(b) A = -|ra 2 , A = - \m 2 (V is equivalent to V5), 

(c) A = — |m 2 , A = — \rn 2 (V is equivalent to Vq). 

2. Three simple Darboux points, when A = 0, and A(m 2 + A) 7^ 0. There are no integrable 
families here as 1^3 = 0. 

4 o 

3. Two simple Darboux points, when either A = and A(m + A) 7^ 0, or A = A = 0. 
Again, no integrable families are present here because Z4 2 = 0. 

4. A triple Darboux point, when A = — m 2 . Additionally there is a simple Darboux point 
when A 7^ 0. The potential is equivalent to V3 and is only integrable when A = —m 2 . 

There are two immediate implications that follow. Firstly, that the main system itself with 
k = is only integrable in those four cases, and the respective first integrals are known, as 
given in the table. Secondly, as was shown in [5] those cases are the only ones which could be 
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integrable when k 7^ 0. This happens because the integrability of the full potential implies the 
integrability of the homogeneous parts of the maximal and minimal degree (the latter is trivially 
solvable in our case). 

As the table shows, when the potential is equivalent to V3 (or, to be precise, its integrable 
subcase) or V4, the second first integral is known; but V5 and V§ only have known integrals with 
zero curvature. And as was shown in [2], for k = 1, the values of A and A are those of V5 or 
Vq forbid integrability. This is easily extended to the k = — 1 case, since after the change of 
variables 

Pj ^e-^ 4 Pj , j = 1,2, (16) 

we obtain a system with the sign of k changed, but the ratios m 2 /A and m 2 /X the same. Thus, 
concerning the conjecture of the quoted paper, our results for k 7^ enable us to state, that it 
is true, when rational integrability is considered. 

However, the above considerations assume that the energy value is generic, so that the 
particular solution is a non-degenerate elliptic function. As stressed in the first part, this does 
not preclude the existence of an additional first integral on the physically crucial zero-energy 
level. 



5 Integrability on the zero-energy hypersurface 

We choose not to investigate the Darboux points, but the variational equations directly, as 
they are considerably simpler in this case. The Hamiltonian equations of ([6]) are 



(17) 



91 = Pi, Pl = -kqi + m 2 qiql - Aqf, 
q 2=P2 , p 2 = -kq 2 + m 2 q 2 q 2 -\ql), 

and they admit three invariant planes as was shown in [9]. They are 

Ilfc = {(qi,q2,Pi,P2) e C 4 I q k = Ap k = 0}, k = 1,2, 

n ^4, i 2 m 2 + A ( 18 ) 

II3 = {(.qi,q2,Pi,P2) G C \q2 = ctqi Ap 2 = -api\, a = 

m z + A 

Obviously two particular solutions are 

{qi =Pi=0, q 2 = q 2 {v), P2 = q^iv)}, ® = \(pI + k ql + > 
fe = P2 = 0, qi = 31(77), pi = (/i(rf)}, = i (+p\ + kq\ + ^qf^j , 



(19) 



and in order to find the third particular solution we make a canonical change of variables 

(qi,q2,Pi,P2) T = B(Q U Q 2 , P u P 2 f (20) 
where symplectic matrix B has the block structure 
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and a and b are defined by 



m 2 + A , / m 2 + A 



2m 2 + A + A' V2m 2 + A + A' 



(22) 



Let us introduce five quantities 



(23) 



(25) 



ai=2m 2 + A + A, a 2 = 3AA + 2m 2 (A + A) + m 4 , a 3 = y/(X + m 2 )(A + m 2 ), 
a 4 = A 2 + A 2 - AA - m 4 , a 5 = AA - m 4 . 

Then, in the new variables, Hamiltonian ([5]) has the form 

H= X - [P 2 + P 2 + k(Q\ + Q 2 )] + J- [a 5 g 4 + 2a 2 Q 2 Q 2 . + 4(A - A)a 3 QiQ2 + «4<9|] • (24) 

and the Hamiltonian equations read 

Qi = Pi, A = -A;Qi - — [a b Q\ + aaQiQ 2 . + (A - A)a 3 Ql] , 

Q2 = P 2 , P 2 = -feQ 2 - — [a 2 Q\Q2 + 3(A - A)a 3 QiQ| + CK4Q2] • 

ai 

Thus, the third particular solution can be seen to be 

{Q 2 = P 2 = 0, Qi = Qi(v), Pi = +Q'i (v)h = l( p ? + kQ * + ^t 1 Q/ t) ■ ( 26 ) 

Of course, this is only valid for a\ ^ 0. We investigate what happens when A + A = —2m 2 at 
the end of this section. 

Normal variational equations along those three solutions (in the position variables) are 

£"(v)= [-k + m 2 q(r,) 2 ] £(77), 

C"( V ) = [-k + m 2 q( V ) 2 ]^r,), (2?) 



-k - —q{r}) 2 



where q{rj) is one of {qi(r]), (72(77), Qi(v)}i depending on the respective particular solution. 

We will consider the k = case first. Changing the independent variable to z = q(f]) 2 , all 
the NVE's are reduced to the following 

2z 2 f(z) + 3zti'(z) -\i£(z) = 0, (28) 

whose solution is 

£(*) = z -d±Vl+8A-)/4 ) (29) 

where we have introduced three important quantities 

Al =-^, A 2 = -f, A 3 = ^. (30) 
A A 05 
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Note, that if any of A, A or 05 is zero, the corresponding particular solution is constant and 
cannot be used to restrict the problem's integrability. Thus, we are left with the E = k = case 
as potentially integrable. 

When we assume k 7^ 0, or equivalently k 2 = 1, and introduce the same independent variable 
z as before, the NVE's read 

2z 2 (Az + 2k)C"(z) + z(3Az + 4Je)g(z) + (m 2 z - k)((z) = 0, 

2z 2 (Xz + 2k)£"(z) + z(2,\z + 4h)g(z) + (m 2 z - k)£(z) = 0, ^ 

2z 2 (^z + 2k\ g'(z) + z U^-z + 4k\ ?(z) - (^z + kj £{z) = 0. 

First, let us observe that unlike in the previous case, when any of A, A or a§ is zero, the system 
is not integrable. This happens, because the NVE's then becomes the Bessel equation 

sH"(s) + S ^(s) + (s 2 -n 2 )as) = 0, (32) 



with n = 1 and in a new variable s = my/ ' z/k (for the first two) or s = my/ ' —2z/k (for the 
third). Such equation is known not to posses Liouvillian solutions [8]. Together with the results 
of the previous section this leads us to the following 

Lemma 1 System )[3j) considered on the zero or generic energy hypersurface with k 2 = 1 is 
not integrable when A or X is zero. Additionally for A + A / —2m 2 , it is not integrable when 
AA = m 4 . 

Assuming that none of those constants is zero, we rescale the variable z in the three equations 
with 

2k 2k 2kai , , 

z^ — -z, z^ — -z, z^ z, (33) 

A A 05 

respectively, so that all three are transformed into a Riemann P equation of the form 

«(*)=<>, (34) 



«•<,) + (iziz£ + iz!z!] 4 . w + 



s#_ 7 y pp 1 - 55' - 7 y 

z 2 + (z- l) 2 + z(z-l) 
with the following pairs of exponents (5,5'), (7,7')) (ft, ft') at its singular points 

(\ + \^ 1 + 8X ^\~\^ 1 + 8X ^ i = l,2,3. (35) 

Using Kimura's results on solvability of the Riemann P equation [7j we check when the difference 
of the exponents give us integrable cases, and find that the parameters must belong to the 
following families 

A . = ^+i) ) j. gZj i=12 3 (36) 

These polynomials in ij are invariant with respect to the change — 1, so it is enough to 

consider non-negative values only. Furthermore, Ai and A2 cannot be equal to zero, as m 2 7^ 0, 
so l\ and I2 need to be strictly positive. 
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This result can be refined still. First, let us notice, that Aj are not independent, for a, are 
functions of A and A. We find the relation between them to be 



with the exception of Aj = 1. Ai and A2 cannot both be equal to 1, as that would mean = 0. 
If only one of them, say Ai is 1, then necessarily A3 = 1 and the only possibly integrable cases 
are those with A2 satisfying (|37[) with I2 > 2. The same holds when Ai and A2 are interchanged. 

When l\ and I2 are taken to be grater than 1, Ai and A2 are positive, so the relation (J3T1) 
requires that 2/{\% — 1) is negative. This only happens for I3 = and it follows that h = h = 2, 
which is exactly the first known integrable case. Since l/(Ai — 1) and 1/(A2 — l)are positive and 
tend to zero monotonically as > 2 tends to infinity, there are no other solutions, and no other 
integrable sets of parameters. 

Finally, we turn to see what happens when A + A = —2m 2 . This is equivalent to 

11 

- + - = 2, (38) 
Ai A2 

and the same two conditions of (|36p because the first two variational equations can still be used. 
It is straightforward to check that the only integer solution of 

1 1 

h(h + l) + h{h + l) ~ 1 (39) 
is l\ = I2 = 1, which we recognise as the second case of our table. 



6 Conclusions 

Bringing the results of both parts together we can state the following properties of the 
system. For the conformally coupled scalar fields: 

Theorem 1 System with a generic energy hypersurface is integrable if, and only if, 

1. k = 0, and its parameters belong to the four families listed in the table. Otherwise there 
exist no additional, meromorphic integral. 

2. k 2 = 1, and its parameters belong to the first two families of the table. Other than that, 
there exist no additional, rational first integrals. 

The second part of the above theorem can be strengthened to meromorphic integrals, although 
not for all values of the parameters, as described in [2|. 

Theorem 2 // system when restricted to the zero energy hypersurface, is integrable, then 
either 

1. k = 0, or 
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2. k 2 = 1 and its parameters belong to the first two families of the table, or 

3. k 2 = 1 and one of {\\, A2} is equal to 1, and the other satisfies the condition $36\) with 
k > 2. 

Otherwise, the system is not meromorphically integrable. In particular this means, that for 
k 2 = 1 if at least one of A and A is zero, the system is non-integrable. 

For the minimally coupled scalar fields, given in [12] by system (14): 

Theorem 3 For A = 0, if the system is integrable then necessarily E = k = 0. 

Theorem 4 When A / 0, if the system is integrable on a generic energy level then either 

1. k = and 9 - 4m 2 /L = (2n + I) 2 , or 

2. 9 — 4m 2 /L = (2n) 2 (regardless ofk), 
for some n £ 7L. 

Theorem 5 For the zero energy hypersurface, and provided that A / 0, if the system is inte- 
grable then either 

1. k = 0, or 

2. 9- Am 2 /L = (2n + l) 2 , n G Z. 

Of course, depending on the properties of the first integrals, we might get quite different re- 
sults, and the requirement of meromorphicity or rationality is still very restricting. As described 
in the introduction, this leaves open the question of existence of real-analytic first integrals. 
Also we recall that physically the scale factor a cannot even assume negative values, and some 
authors argue that when cosmological (instead of conformal) time is used, the evolution is not, 
in essence, chaotic [3]. Thus, we would like to stress that Liouvillian integrability is a mathemat- 
ical property of the system, and often the methods used to study it require the complexification 
of variables. This means that when restricted to the narrower, physical domain, the dynamics 
might be much simpler. And in particular we might be interested in a particular trajectory 
whose behaviour is far from generic. It is no surprise then, that the dynamics of our system 
when restricted to a > might appear regular. It should still be noted that the notion of chaos, 
although frequently associated with integrability, has not yet been successfully conflated with 
it. And that regular evolution is not necessarily integrable. 
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Appendix A. Massless field 

For m = we can separately solve for each variable, so that we have 



E 1 = ~a 2 - ^ka 2 + ^Aa 4 , 
E 2 = \^ + \f 2+ \k, 2 + \x,\ 



(40) 



with Ei + E 2 = E being the total energy. The first of these is immediately solved, when we 
substitute v\ = a 2 to get 

vf = 2Avf - Akvf - 8E lVl , (41) 

whose solution is 

2 2k 

vi(v) = -£p(v-m;92,g3) + 2£, (42) 

with 771 the integration constant and 

g 2 = ^k 2 + 4AE 1 , g3 = ^-k 3 + ^kAE 1 . (43) 

Of course, when A = the Weierstrass function p reduces to a trigonometric function. 
Similarly, for the other variable, we substitute v 2 = 4> 2 and obtain 

v\ = -2Xv 3 -4kv 2 + 8E 2 v-4lo 2 , (44) 

whose solution is 

2 2k 

V2(v) = -jp(v- m; 54,55) - g^, (45) 

where 

54 = -k 2 + AXE 2 , g 5 = ^k 3 + -kXE 2 + X 2 iv 2 , (46) 

and 772 is the integration constant. As before, for A = the solution degenerates to trigonometric 
functions. 
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